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SUBMILLIMETER NARROW EMISSION LINES FROM 
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ABSTRACT 

A spectral-line survey of IRC-l-10216 in the 345 GHz band has been undertaken with the Submil- 
limeter Array. Although not yet completed, it has already yielded a fairly large sample of narrow 
molecular emission lines with line-widths indicating expansion velocities of ~4 km s~^, less than 3 
times the well-known value of the terminal expansion velocity (14.5 km s~^) of the outer envelope. 
Five of these narrow lines have now been identified as rotational transitions in vibrationally excited 
states of previously detected molecules: the v=l, J=17-16 and J=19-18 lines of Si'^'^S and ^^SiS and 
the v=2, 3—7-6 line of CS. Maps of these lines show that the emission is confined to a region within 
~60 AU of the star, indicating that the narrow-line emission is probing the region of dust-formation 
where the stellar wind is still being accelerated. 

Subject headings: stars: individual (IRC-t-10216) — stars: late-type — circumstellar matter — sub- 
millimeter — radio lines: stars 



1. INTRODUCTION 

Massive circumstellar envelopes of Asymptotic Giant 
Branch (AGB) stars are believed to be a major contrib- 
utor of molecules and grains to the interstellar medium. 
IRC-H10216 (CW Leo), at a distance of 150 pc, is the 
archetypal AGB carbon star with a high mass-loss rate 
■ (> 10-^ M0 yr-i) (Young et al. 1993; Crosas & Menten 
' 1997). Owing to its proximity, this star is an ideal target 
for detailed studies of physical and chemical processes in 
AGB circumstellar envelopes (e.g. Olofsson et al. 1982). 
Nearly 60 molecular species have been discovered in its 
circumstellar shell from single-dish spectral-line surveys 
(Kawaguchi et al. 1995; Avery et al. 1992; Groesbeck et 
; al. 1994; Cernicharo et al. 2000; He et al. 2008). 

Interferometric mapping allows us to check predic- 
tions of abundances of various molecular species as a 
function of radius in the circumstellar envelope. "Par- 
ent" molecules, such as CO, C2H2, CS, HCN, & SiS, 
are formed in the stellar atmosphere in thermo-chemical 
\ equihbrium (Tsuji 1964, 1973). Once they get levitated 
to a distance from the star at which the density is too 
' low for chemical reactions, their abundances "freeze out" 
at the values prevailing in that region (McCabe et al. 
1979). At around that distance (~ 20 AU), the tem- 
perature has dropped below the dust condensation value 
(~ 1200 K) and dust grains start forming (Monnier et 
al. 2000). Radiation pressure accelerates the grains and, 
by friction, the molecules too, in an outflow which, for 
IRC-l-10216, reaches a terminal velocity of 14.5 km s^^. 
In the outer parts of the expanding envelope, a rich 
carbon-dominated photo chemistry driven by the ambi- 
ent ultraviolet field produces species such as CN and C4H 
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whose mostly optically thin emission can be observed as 
ring-like distributions with radii, at a few times 10^^-10^^ 
cm depending on the chemical reactions at work (see the 
reviews of Glassgold 1996; Ziurys 2006). 

At submillimeter wavelengths, we can observe lines 
requiring elevated excitation conditions, allowing us to 
probe the physical conditions in the inner circumstellar 
envelope (radius < 10^^ cm) where densities and tem- 
peratures are relatively high (temperatures ^100-1000 
K, column densities ^ 10^^-10 cm^^) (Keady & Ridg- 
way 1993). 

We have begun a spectral-line survey of IRC+10216 
with the Submillimeter Array^ (SMA, Ho et al. 2004). 
Previous 345 GHz single-dish line surveys include those 
done with the JCMT (Avery et al. 1992) and CSO 
(Groesbeck et al. 1994) in the frequency ranges of 339.6- 
364.6 GHz and 330.2-358.1 GHz, respectively. Our SMA 
line survey will eventually cover the frequency range of 
300-355 GHz with higher sensitivity and spatial reso- 
lution than previous surveys; about 40% of the survey 
has so far been completed. The frequency range 300-330 
GHz will be observed for the first time. Here we present 
several new results, including the discovery of many lines 
having narrow widths, and detections of vibrationally ex- 
cited rotational transitions in SiS, ^^SiS, Si^'^S and CS. 
A full account of the observed lines will be presented on 
completion of the survey. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed IRC-H10216 with the SMA in 2007 Febru- 
ary in the subcompact configuration with baselines from 
9.5 m to 69.1 m in the frequency range of 300 to 355 
GHz. To follow up on some of the detected narrow lines 
with higher angular and spectral resolution, we repeated 
the observations at 337.5 GHz with the SMA in the ex- 
tended configuration on 2008 February 19. The base- 

^ The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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line lengths in this configuration range from 44.2 m to 
225.9 m. The synthesized beam sizes were 3" x 2" and 
0."8 X 0."6 in the subcompact and extended array obser- 
vations, respectively. Table 1 summarizes the observa- 
tional parameters on the five tracks of observations which 
are relevant for the data presented here. The duration of 
each track was from 7 to 9 hours. The phase center was 
a(2000) = 09''47'"57.38% (5(2000) = -hl3°16'43."70 for 
all observations. All the tracks in subcompact configura- 
tion were carried out in mosaiced mode, with 5 pointings 
with offsets in RA and DEC: (0", 0") and (±12", ±12"). 
The extended configuration observations were carried out 
with single pointing toward IRC+10216. Titan and the 
quasars 0851+202 and 1055+018 were observed every 20 
minutes for gain calibration. The spectral band-pass was 
calibrated using observations of Mars and Jupiter. Ab- 
solute fiux calibration was determined from observations 
of Titan and Ganymede. 

The visibility data were calibrated using the Miriad 
package (Sault, Teuben & Wright 1995). The mosaiced 
images are deconvolved using the Miriad task mossdi; 
the resulting synthesized beams are summarized in Ta- 
ble 1. Maps of continuum emission show the peak to 
have a position offset of (Aa, AS) « (0."7, 0."2) from 
the phase center position quoted above. The absolute 
position measurements in the continuum emission are es- 
timated to be accurate to ^ 0."1. Taking into account 
the proper motion of IRC+10216 of (d, S) « (26, 4) mas 
yr~^ determined by Menten et al. (2006) our position is, 
within the mutual uncertainties consistent with that de- 
termined by those authors. The continuum emission was 
unresolved at the highest angular resolution of ^ 0."8. 
The integrated continuum flux density was 0.84 Jy at 
301.1 GHz and 1.17 Jy at 337.5 GHz, with an uncertainty 
of about 15% in the absolute fiux calibration. All the 
spectra shown below were produced by integrating the 
continuum-subtracted line intensity in a 2" x 2" rectan- 
gle centered on continuum peak (by means of the Miriad 
task imspec). 

3. EXPANSION VELOCITIES 

Within a 3" beam a total of 92 lines were detected 
in the first phase of the SMA line survey of IRC+10216 
at the central position. Truncated parabolic line profiles 
were fitted using the CLASS package (using the shell 
model for the line-profile) with one of the fitted param- 
eters being Vexp, the expansion velocity of the circum- 
stellar shell. From this sample, 25 lines have Vexp < 
7 km s^^and of these, 12 are as yet unidentified. Sev- 
eral are tentatively identified to be lines of salts such as 
KCN, NaCl and NaCN. Some of the lines may result from 
known molecules in vibrationally excited states, such as 
the V3=l 157^9 - 147,8, 157,8 - 147,7 doublet of SiCC at 
345727.3 MHz. Figure 2 shows a sample spectrum to- 
ward IRC+10216 observed on 2008 February 9 over a 
2 GHz wide band centered at 337.5 GHz. A compari- 
son with the line survey of Groesbeck et al. (1994) (see 
their Figure 1) shows that only the C'^^S J=7-6 fine at 
337.396 GHz was detected in their observations. All of 
the new narrow lines were missed in this previous survey 
due to poorer sensitivity. Here we present results on the 
lines which are securely identified as rotational transi- 
tions in vibrationally excited states of molecules that are 
well known to be abundant in IRC+10216's envelope. 
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Fig. 1. — Distribution of expansion velocities derived from lines 
detected in the circumstellar envelope of IRC+10216. The uncer- 
tainty in Vexp is ^0.2 km s~^(1(t). White bars with bold outlines 
are from the recent line survey of He et al. (2008) which is repre- 
sentative of all single-dish line surveys toward IRC+10216 at cm 
to submm wavelengths. The bin at 14 km s~^is shown truncated 
here; it actually consists of 170 lines. Grey bars represent SMA ob- 
servations, showing a new population of narrow lines which peaks 
at ~4 km s~^. 

In previous single-dish line surveys of IRC+10216 (Cer- 
nicharo et al. 2000; Kawaguchi et al. 1995; Avery et al. 
1992; Groesbeck et al. 1994) Vexp is not tabulated, but 
it can be seen from published spectra to be ^^14. 5 km 
for all lines. From the latest published line survey 
by He et al. (2008) (see their Table 9.), we can plot a 
distribution of V^xp which is shown in Figure 1 as empty 
bars with bold outlines. This histogram peaks at 14 km 
s~^. The distribution of Vexp from our line survey is 
shown in Figure 1 as grey bars. Our line survey shows 
a peak in the same bin of 14 km s~^but reveals a signif- 
icant number of narrow lines with velocities around '^4 
km s~^. Both histograms show a continuous distribution 
of expansion velocities between these two peaks at 4 and 
14 km s~^ 



Lines with Vexp < 10 km s~^ from IRC+10216 have 
been reported by Highberger et al. (2000) and were as- 
signed as vibrationally excited SiS and CS lines. He et 
al. (2008) found four lines with Vexp=7~10.2 km s~^ (see 
their Table 14) - all from vibrationally excited SiS. Nar- 
row maser lines in IRC+10216 from SiS and HCN were 
reported by Fonfrfa-Exposito et al. (2006) and Schilke 
& Menten (2003), respectively Savik-Ford et al. (2004) 
detected OH - which has a narrow width (5.8 km s~^), 
but this line appears at —37 km s~^, blue-shifted with 
respect to the systemic velocity of the star of —26 km 
s^^. All the narrow lines we have observed are centered 
at the systemic velocity of about -26 km s^^to within 
~2 km s~^. Infrared molecular line profiles in the 10 /im 
band were observed and analyzed by Keady & Ridgway 
(1993). They proposed (see their Figure 3) an expansion 
velocity as a function of radius with: (1) Yexp = 3 ~ 4 
km s^^ from l'^ 8 R,, (2) rapid acceleration from 4 to 11 
km s^^ from 8 to 10 R* and (3) slower acceleration from 
11 km s^^to the terminal velocity of 14 km s~^ from 10 
to 20 R* . The histogram of expansion velocities (Figure 
1) peaks at ~4 km s~^ and ~14 km s^^, consistent with 
this velocity structure. 
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Fig. 2. — A sample spectrum toward IRC+10216 showing several 
examples of narrow lines. Over this frequency range, the only line 
detected in a previous line survey (Grocsbcck et al. 1994) was the 
C^^S J=7-6 line at 337.396 GHz. This line has a Ve^p =13.9 km 
s~ ^ . All the other lines are new detections and several of them are 
narrow (Vexp < 7 km s~^). 



4. VIBRATIONALLY EXCITED SIS, ^SSIS AND SI^^S 

SiS is an abundant molecule in IRC+10216, one of the 
parent molecules foimd close to the star (Glassgold 1996), 
and an important source of Si for the formation of silicon 
carbides such as SiC, SiC2 and SiC4 (e.g., McCarthy et 
al. 2003). Figure 3 (top) shows the spectrum of SiS v=l 
J=19-18 emission at 343100.98 MHz (rest frequency), ob- 
served on 2007 February 12. This line was previously de- 
tected (but not mapped) in the single-dish line-survey of 
Groesbeck et al. (1994), where it is significantly weaker 
than here because of poorer sensitivity (we have longer 
integration time as well as larger collecting area with the 
SMA) . Figure 3 (bottom) shows a map of the integrated 
intensity emission, over the velocity interval of -40 to -10 
km s^^. The total integrated flux density is 69.1 Jy km 
s~^. The deconvolved source size is l.l"x0.5" with P.A. 
of —27.6°. Assuming a size of 0."5 x 0."5, we estimate a 
lower limit to the brightness temperature of 200 K. 

The emission appears to have an extended weak feature 
towards the southeast. High angular resolution near-IR 
adaptive-optics images of IRC+10216 show evidence of 
azimuthally asymmetric structures in the inner circum- 
stellar envelope, over angular scales of ^ 2" (Menut et 
al. 2007). The 3" angular resolution of our observations 
is insufficient to allow a detailed comparison between the 
submillimeter line emission maps and near-IR images. 

The expansion velocity of the SiS v=l J=19-18 line is 
10.6 km s~^. This is an example of a line of intermediate 
expansion velocity (see Figure 1), in the rapidly accel- 
erating zone of the envelope where presumably dust has 
already formed. 

Figure 4 shows spectra of Si^^S and ^^SiS v=l, J=17- 
16 and 19-18 lines, detected toward IRC+10216 for the 
first time. This emission is unresolved and the source- 
size upper limits are shown in Table 2. Assuming the size 
of the emitting region is 0."2, we estimate a brightness 
temperature to be 140~200 K, which should be consid- 
ered as lower limits. We estimate a column density of 
7x10^'' cm^^ and abundance of 7x10^^ for ^^SiS for an 
assumed excitation temperature of 550 K. 




Velocity (VIsr km/s) 




Fig. 3. — Upper: Spectrum of SiS v=l J=19-18 line obtained 
from a 2" square centered on the star, with Veip=10.6 km s^^. 
For comparison, the SiS, v=0 ,1=19-18 line is shown by the dashed 
line, (scaled down by a factor of 10). The v=0 line has an ex- 
pansion velocity of 13.5 km s^^. Note that the intensity val- 
ues for this line arc scaled down by a factor of 10 in this fig- 
ure. Lower: Integrated intensity emission from the SiS v=l J=19- 
18 line. The contours levels arc -5, 5, 10, 20, 40, 80, 120 X 
0.45 Jy/beam km/s. The coordinate offsets arc with respect to 
q(2000) = 09'' 47™•57.43^ 5(2000) = +13°16'43."98.Tlic synthe- 
sized beam is shown in the lower right corner. See also Table 2. 

From the observed integrated intensities of SiS, ^^SiS 
and Si^^S v=l, J=19-18 emission, we derive the iso- 
topic abundance ratios p^Si/^^Si] = 15.1 ± 0.7 and 
pg/34g] ^ 19.6±1.3 (uncertainties are 1<t). For compar- 
ison, previously published values from single-dish obser- 
vations (Kahane et al. 1988; He et al. 2008) are 20.2 ±2 
for p^S/^^S](in good agreement with our estimate here) 
and 18.7 ± 1 for [Si/^^Si] (marginally larger than our 
value). A possible cause for disagreement might be non- 
negligible optical depths. 

5. VIBRATIONALLY EXCITED CS 

Vibrationally excited CS radio emission in IRC+10216 
was first detected by Turner (1987) in the v=l, 3=2 1 
and 5 4 transitions. These lines were re-observed along 
with new detections of the J =3-2, 6-5 and 7-6 tran- 
sitions by Highberger et al. (2000). The J=7-6 emis- 
sion reported here is the first detection in the v=2 state 
(Figure 5). The triangular line-profile is indicative of 
spatially unresolved emission from accelerating gas (Bu- 
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Fig. 4. — Spectra of vibrationally excited rotational lines in iso- 
topes of SiS obtained from the imaged data cube. The flux density 
was calculated over the area enclosed in a 2 " square centered on 
the star. See also Table 2. 
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Fig. 5.— Spectrum of CS v=2 J=7-6 emission at 337912.19 GHz. 
See Figure 4 caption and also Table 2. 

jarrabal et al. 1986). In subcompact configuration ob- 
servations with angular resolution of ~ 3", the emission 
appears highly concentrated and unresolved. There is no 
emission detected at the angular radius of ~ 12" (where 
several other species have shown a peak in abundance in 
previous interferometric maps). The observations of this 
line were repeated in the extended configuration of the 
SMA at a beam size of 0."8, confirming that the emis- 
sion is unresolved. The deconvolved source size is < 0."2. 
The lower limit for brightness temperature is 237 K. As- 
suming an excitation temperature of 550 K, we estimate 
a column density of 7 x 10^"^ cm~^ and a lower limit for 
CS abundance with respect to H2 of 9.3 x 10^^ within a 
radius of - 4.5 x lO" cm 7i?,). 



Previous interferometric observations of the v==l, J=5- 
4 line show the emission to be within a radius of 0."35 
(~10 R, or ~ 7 x 10" cm) (Lucas and Guelin 1999). 
Young et al. (2004) derived a lower limit on the CS 
abundance to be 3.4 x 10^^ within ~ 34i?,. From single- 
dish observations of CS v=l, J=3-2, 6-5 and 7-6 lines, 
Highberger et al. (2000) estimate an abundance of 3- 
7xl0~^ relative to H2. From comparison with the CS 
radial abundance predicted by chemical models, we con- 
clude, as did Young et al. (2004), that the model of 
Millar et al. (2001) predicts too low a value (-- 10~") 
at the radius of ^ 10^^ cm. In this model, assuming 
CS to be a parent molecule (see Figure 1, right-panel, in 
Millar et al. 2001), the initial abundance is still lower 
by about an order of magnitude relative to that derived 
from submillimeter observations. Moreover, the drop in 
CS abundance with radius (owing to the production of 
other sulphur-bearing molecules), is too small, and not 
consistent with the compact distribution of CS seen in 
the SMA observations. There is better agreement in a 
more recent study of non-equilibrium chemistry of the 
inner wind, which takes into account shocks induced by 
stellar pulsation (Cherchneff 2006), although these mod- 
els seem more relevant for S stars (C/0~ 1). 

The CS v=2, J=7-6 transition requires extreme ex- 
citation conditions since it corresponds to an energy 
E„/k=3707 K. As noted by Highberger et al. (2001), 
even for the v=l line, coUisional excitation with II2 would 
require very high gas densities (^ 1 — 5 x 10^"* cm~'^. This 
line emission is most plausibly excited by 8 /im stellar 
thermal radiation. 

6. CONCLUSIONS 

Preliminary results from the SMA line survey of 
IRC-l-10216 have yielded a population of narrow lines 
with expansion velocities of ^ 4 km s^^. About half of 
these can be assigned to vibrationally excited rotational 
transitions of abundant species such as CS, SiS and their 
isotopomers. The emission is found to occur in a very 
compact region smaller than 0."2 around the star. This 
is thought to be the region where dust is forming in the 
envelope and in which the material has just begun ac- 
celerating and has yet to attain the terminal velocity of 
~14 km s"^ The CS v=2, J=7-6 line is most hkely 
radiatively excited, because coUisional excitation would 
require an unrealistically high gas density and abundance 
of CS. 

It is a pleasure to thank Ray Blundell for his 
help and support present on the SMA IRC+10216 
line-survey. We thank Mark Gurwell, Thushara 
Pillai and Jun-Hui Zhao for helpful discussions on 
SMA data reduction. This research has benefitted 
from the Cologne Molecular Spectroscopy Database 
(Miiller et al. 2001; Miiller et al. 2005) 
( http: / /www. astro.uni-koeln.de /site / vorhersagen/) , and 
the ALMA group's spectral line catalog website: 
] http : / / www .splatalogue.net (Remijan et al. 2007) 
and the CASSIS (Centre d'Analyse Scientifique de 
Spectres Infrarouges et Submillimetriques) software 
( |http://cassis.cesr.£r| . 
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TABLE 1 
Summary of observations 



Date 


SMA 


Tuning 


Synthesized 




T225GHZ 






configuration 




beam 






(SSB, K) 


2007 February 7 


subcompact 


299.1 GHz (LSB) 


3. "2 X 2. "4, P.A.= 


-4° 


0.08 


180-270 


2007 Fcburary 8 


subcompact 


301.1 GHz (LSB) 


3."3 X 2. "4, P.A.= 


-6° 


0.09 


160-230 


2007 Fcburary 9 


subcompact 


337.5 GHz (LSB) 


3."0 X 2. "4, P.A.= 


-8° 


0.08 


180-350 


2007 Fcburary 12 


subcompact 


334.4 GHz (LSB) 


3."0 X 2. "2, P.A.= 


-3° 


0.05 


130-260 


2008 Feburary 19 


extended 


337.5 GHz (LSB) 


0."8 X 0."6, P.A.=- 


88° 


0.03 


100-250 



TABLE 2 
Summary of narrow lines 



Species 


Transition 


Rest frequency^ 


Catalog frcq^. 


^ exp 


Peak 


Integrated 


Deconvolved 














flux density^ 


flux density 


size 








(MHz) 


(MHz) 


(km s^-') 


(Jy) 


(Jy km s~^) 




(K) 


OS 


v=2 J=7-6 


337913.246 ± 0.18 


337912.189 


5.1±0.1 


0.78 


7.0 


< 0."2 


>237.0 




v=l J=17-16 


298630.441 ± 0.69 


298629.989 


4.5±0.4 


0.36 


2.6 


< l."9 X 0."2 


>138.7 


Si34S 


v=l J=19-18 


333733.581 ± 0.32 


333731.998 


5.7±0.1 


0.49 


4.5 


< 2. "3 X 0."5 


>183.3 


29SiS 


v=l J=17-16 


301390.406 ± 0.25 


301388.939 


4.8±0.2 


0.52 


3.8 


< l."9 X 0."8 


>198.3 


29SiS 


v=l J=19-18 


336819.842 ± 0.32 


336814.954 


7.5±0.1 


0.50 


4.6 


< 0."9 X 0."6 


>152.9 



^ Assuming a systemic velocity of —26.2 km s ^.^ Prom the Cologne Database of Molecular Spectroscopy.^ Typical uncertainty is 0.1 
Jy.^ Lower limit in brightness temperature assuming source size of 0."2 x 0."2 
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